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ABSTRACT 

Deep drawing is a complex process that is influenced by mrious parameters like tooling con/igurations, blank 
con jigurations, material properties, and f orming conditions. Limit Drawing Ratio and limit strains are two parameters 
to indicate the draw ability of the material. In this paper, the effect of various parameters like die and punch corner 
radii, clearance, blank holding force, friction coefficient between die and punch, punch diameter on limit drawing ratio 
is investigated andpresented by using an explicit finite element code LSDYNA 3D. The model is validated by comparing 
force obtained by the existing experimental setup with the one obtained with simulation. Initial estimation of the effect 
of the parameters listed above and their interactions is found out by using two levels L8 orthogonal array. After 
identifying the most effective parameters, their effect on the LDR is studied andpresented. 
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INTRODU CTION 

Deep drawing process is used to manufacture complicated parts from sheet metal and is used in many 
industries such as automobile, aerospace, appliance and so on. The process design of the deep drawing process 
involves so many parameters such as tooling configurations, blank configurations, material properties, and forming 
conditions, thereby involving a lot of trial and error runs which increase not only the cost of development but also 
the time to market. Therefore, systematical design and process simulations have presented themselves as a vital 
part in today’s competitive systems. 

Drawability is defined as the ability to draw the blanks into cups of desired shape without any defects and 
is difficult to be quantified. One of the parameters which indicate the drawability is the limit drawing ratio, which 
is the ratio of the maximum blank diameter to the cup diameter at the onset of tearing. The limit drawing ratio is 
usually limited by the maximum allowable drawing force at the beginning of the punch nose rounding, when 
punch nose radius is not too great and appreciable stretching over the punch nose done not take place [1] 
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Finite Element Model & Experimental Yalidation 



True Strain 


Figurel: Finite element model Figure 2: True Stress-Strain Diagrams of Specimens 

at 0°, 450 and 900 to Rolling Direction 

Initially an FEA Model (Figure 1) is created in the existing setup of punch, die, blank and blank holder for which 
the comparison of force obtained from the simulation with the experimental result. The existing set up is designed to 
produce a 25mm cup from a 45 mm blank. Based on the symmetry boundary condition, a quarter of the geometry is 
modeled by using a Hyper mesh with LSDYNA template. The punch presses the sheet into the die cavity. In sheet metal 
forming, generally, membrane element or continuum element or shell elements are employed [2]. Since membrane 
elements lack the bending stiffness and the continuum element takes higher computation time, the blank is modeled with 
shell elements of the mid plane with Belytschko formulation and with five integration points. Punch, die and blank holder 
are taken as rigid materials and for obtaining material model of the blank material 


Tensile test specimens of ASTM Standard size are prepared from the commercial pure Aluminum, AA1100, 
sheet. The pieces are cut in the rolling direction, 45° to rolling direction and in transverse to the rolling direction The pieces 
are tested at AMTL (Aero Space Material testing Laboratory), A division of MIDHANI, DRDO on a universal testing 
machine (INSTRAN 4507 MODEL). The true stress and strains are computed up to the ultimate point and presented in 
Figure 2. 



The material properties obtained for the specimens in three directions are presented in table 1. 


Table 1: Material Properties Obtained from the Tensile Test 


s. 

No 

Property 

Yalue 

0 

45 

90 

Average 

1 . 

Yield strength (0.2% Proof) in Mpa 

142.7 

138.9 

134.1 

138.6 

2. 

Ultimate Tensile Strength (Mpa) 

161.3 

154.7 

156.6 

157.7 

3. 

Modulus of Elasticity (Mpa) 

57109 

50908 

54312 

54110 

4. 

Yield Strain (from values) 

0.002697 

0.002926 

0.002667 

0.002763 

7. 

Anisotropic Index (R) 

1.12 

1.06 

1.07 


8. 

Planar Anisotropy 


9. 

Normal Anisotropy 

1.078 



Figure 3: Flow curve of AA1100 
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To check the initial anisotropy, the value of the normal anisotropy 4 R ’ is computed and presented in table 1. It is 

observed that the value of R is slightly more than 1. So initial anisotropy is neglected and the blank is treated to be 
isotropic, which simplifies the material characterization. . Considering the above the material model ‘piece wise linear 
isotropic plasticity model is selected for analysis for which flow curve is to be input. Since the strain is additively 
decomposive, the plastic strain ‘s p ’ can be found out by subtracting elastic strain ‘s e ’ from the total strain ‘s’. s e can be 

found out using the relation £ e —G y !E Where a y is the current stress. On obtaining the plastic strain values the curve 

between true stress and plastic strain is plotted and presented in figure 3. The flow curve obtained in the transverse 
direction is taken as input, as it is noted from the table that the average properties are nearer to the properties in transverse 
direction to the rolling direction. 

The mesh convergence is checked by the convergence of the punch load by the time for various mesh sizes as 
done by Jamal Hematian [3]. 


Experimental validation of the finite element model is carried out by comparing the force obtained by the 
experiment with force obtained in simulation. The experimental setup as shown in Figure.4 consists of a 25T hydraulic 
press interfaced to the computer with the load cells through digital force indicator. The blank holding schema obtained 
from the experimental set up shown figure.5 is applied in the simulation. 




Time in Sec 


Figure 4: Experimental Set up Figure 5: BHF Schema Obtained from the Experiment 



Figure 6: Transient Load v/s Time Diagram obtained from the Experiment 

Since it is spring loaded blank holder the schema should be linear. The same is evident from the Figure 5 and the 
force obtained from the experimental setup is shown in Figure. 6 and by simulation in Figure.7. 
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Figure 7: Transient Load v/s Time Diagram obtained from the Simulation 

The maximum force obtained from the experiment is 1184 Kg. and from the simulation, it is found to be 2697.9N 
i.e 275Kg. Since it is quarter model the actual force obtained from simulation is four times the given value of the 
simulation i.e 1100 Kg. The deviation from the experimental result is 7%. The reason for this discrepancy may be due to 

non-isotropic hardening during forming operation. The failure occurs when o eff = o u i t where o eff and o u i t are von-Mises’ 
effective stress and Ultimate Tensile Strength (110 MPa for AAl 100) respectively. 

In this work the effect of die profile, punch profile, clearance, friction coefficient, punch diameter on limit drawing ratio is 
investigated 

SIMULATIONS 

After the validation, the simulations are planned following DOE techniques [4] 5 

Since each simulation takes a large amount of computational time. The screening test for the above five factors is 
planned by considering the factors at two levels. The screening test levels considered are shown table 2. Generally for 
continuous factors the first level may be the lower of the two values being tested and the second will be the higher. 


Table 2: Levels Chosen for Screening Test 


s. 

Factor 

Level 1 

Level 2 

No 

(Min) 

(Max.) 

1 . 

Punch Profile 

2 mm 

5 mm 

2. 

Die Profile 

3 mm 

10 mm 

3. 

Clearance 

7% 

20% 

4. 

Friction Coeff. 

0.015 

0.45 

6. 

Punch diameter 

lOmm 

40 mm 


Table 3: L8-Orthogonal Array after Assignment of Factors 


Trial 

No 

Column Number 

Rp 

(mm) 

Rd 

(mm) 

c 

% 

P 

D P 

(mm) 

1 

2 

3 

1 

0.015 

10 

2 

2 

3 

20 

0.45 

40 

3 

2 

10 

7 

0.45 

40 

4 

2 

10 

20 

0.015 

10 

5 

5 

3 

7 

0.45 

10 

6 

5 

3 

20 

0.015 

40 

7 

5 

10 

7 

0.015 

40 

8 

5 

10 

20 

0.45 

10 


Clearance is generally from 7% to 14%. H.Gharib et.al [5] has used 20% clearance with die radius up to 10 mm 
for the cylindrical cup of 25 mm diameter with punch and die profiles of 5mm and 10 mm respectively. M. Jurkovic et.al 
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[6] have experimentally found out friction coefficients in deep drawing for yarious surface conditions of tooling at different 
degrees of deformation. They vary from 0.015 to 0.18. In case of dry contact, the friction between and A1 and steel is 0.45 
for sliding [7]. The above range is selected for the analysis for the present work. An L8 orthogonal array is chosen which 
has 7 columns and 8 trials. The orthogonal array after factor assignment is shown in Table 3. 

RESULTS & DISCUSSIONS 


As a sample, von-Mise’s stresses for various drawing ratios for trial 1 and are presented in Figures from Figure 8 
to figure 11. 


The limit drawing ratios for various trials are presented in Table 4. The analysis is being carried out by column 
effect method. Summation at level 1(S1) and the summation at level 2 (S2) are computed and the difference (S1-S2) is 
found out. According to the column effect method, higher the difference between the two sums, higher is the effect of that 
parameter. 



Fringe Levels 
1.041e+02 
9.837e+01 ■ 
9.265e+01 | 
8.693e+01 _ 
8.121 e+01 _ 
7.549e+01 _ 
6.977e+01 _ 
6.405e+01 _ 
5.833e+01 
5.261e+01 fj 
4.689e+01 | 



Fringe Levels 
1.115e+02 
1.034e+02 1 
9.539e+01 | 
8.735e+01 _ 
7.931 e+01 _ 
7.127e+01 _ 
6.323e+01 _ 
5.519e+01 _ 
4.71 5e+01 _ 
3.911e+01 ■ 
3.107e+01 ■ 


Figure 8: Effective Stress (Von- Mise’s) 
Distribution with p=2.5 for Trial-1 


Figure 9: Effective Stress (Von- Mise’s) 
Distribution with p=2 for Trial-1 



Fringe Levels 
9.909e+01 
9.482e+01 ■ 
9.056e+01 | 
8.629e+01 _ 
8.203e+01 
7.776e+01 _ 
7.350e+01 _ 
6.923e+01 _ 
6.497e+01 
6.071 e+01 ■ 
5.644e+01 | 



Time = 0.75456 

Contours of Effective 
max ipt. vaiue 
min=15.6573, at 
max=117.246, at 


Fringe Levels 
1.172e+02 
1.071e+02 1 
9.693e+01 g 
8.677e+01 _ 
7.661 e+01 _ 
6.645e+01 _ 
5.629e+01 _ 
4.613e+01 _ 
3.598e+01 _ 
2.582e+01 ■ 
1.566e+01 | 


Figure 10: Effective Stress (Von- Mise’s ) 
Distribution with p=2.2 for trial-1 


Figure 11: Effective Stress (Von- Mise ? s ) 
Distribution with p=2.1 for trial-1 


Table 4: Limit Drawing Ratios for various Trials 


Trial 

No 

Column Number 

p 

Rp 

(mm) 

Rd 

(mm) 

c 

% 

V 

D P 

(mm) 

1 

2 

3 

1 

0.015 

10 

2.1 

2 

2 

3 

20 

0.45 

40 

1.2 

3 

2 

10 

7 

0.45 

40 

1.6 

4 

2 

10 

20 

0.015 

10 

2.5 

5 

5 

3 

7 

0.45 

10 

1.5 

6 

5 

3 

20 

0.015 

40 

1.8 

7 

5 

10 

7 

0.015 

40 

1.4 

8 

5 

10 

20 

0.45 

10 

2.3 

S1 

7.4 

6.6 

6.6 

7.8 

8.4 


S2 

7.0 

7.8 

7.8 

6.6 

6.0 

S1-S2 

0.4 

-0.8 

-0.8 

0.8 

2.4 
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From Table 4, it is clear that the cup diameter has a pronounced effect on LDR and the effect of die profile, 
friction factor and % clearance is considerable. 



Figure 12: Interaction Graph between Punch Diameter and Friction Factor 

Plotting method is used to find out the combined effect of the parameters considered above was taken two at a 
time. For this each pair of parameters at levels 1-1, 1-2, 2-1, and 2-2 are summed up and plotted as shown in Figures 12 to 
16. According to plotting method, higher the difference in slope, higher is the interaction between the two parameters 
chosen. 



Figure 13: Interaction Graph between Punch Figure 14: Interaction Graph between Punch 
Diameter and % Clearance Diameter and Die Radius 




Figure 15: Interaction Graph between Figure 16: Interaction Graph between Friction 
Die Profile and % Clearance Factor and % Clearance 

From the Figure 15, it is evident that there is inter-dependency of die profile and % clearance in effecting the 
LDR. There is a little interaction of punch diameter with die profile and % clearance (Figure 13& Figure 14). Friction 
factor independently effect the LDR (from Figure 12 & Figure 16) 

CONCLUSIONS 

The following conclusions are drawn from this work 

• The most significant parameter which effect the LDR is punch diameter 

• There is a considerable effect of die profile, clearance and friction factor on LDR 
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• Variation of the punch profile will not affect the LDR much 

• There is a strong interaction of die profile and % clearance in affecting the LDR 

• A little interaction is observed for the punch diameter with die profile and clearance. 

• No interactions are found for friction factor. 
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